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Q\ , Abstract: We consider several ways of how one could classify the various types of soliton solutions 
related to nonlinear evolution equations which are solvable by the inverse scattering method. In 
| | '' doing so we make use of the fundamental analytic solutions, the dressing procedure, the reduction 
. technique and other tools characteristic for that method. 
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1. Introduction 

It is our impression that the question in the title has not been answered satisfactorily even for some 
of the best known type of soliton equations such as the iV-wave equations, the multicomponent 
NLS equations and others. 

We use the term 'soliton solution' as a special solution to a given nonlinear evolution equation 
(NLEE) which is solvable by the inverse scattering method 24,4 . That means that the NLEE 
allows Lax representation: 

[L(A),M(A)]=0, 

where L(X) and M(X) are two linear differential operators. In what follows we take them to be 
first order matrix differential operators 

Lip(x,t,X) = id x ip + U(x, t, X)ip{x, i, A), 
Mip(x,t,X) = id t ip + V(x,t,X)ip(x,t,X). 

The one-soliton solutions are related to one or a set of several discrete eigenvalues of the Lax 
operator L. Therefore one first has to study the different configurations of discrete eigenvalues of 
L, see [15]. The next step in classifying the types of one-soliton solutions is related to the study of 
their internal degrees of freedom. 

In order to make the problem not too difficult we will specify L to be the generalized Zakharov- 
Shabat system: 

L(X)i)(x, X) = id x i> + (q(x) - XJ)il>{x, X) = 0. 

where we take the potential q(x,t) to be n x n matrix-valued smooth function of x tending to 
zero sufficiently rapid as x — > ±00. We also restrict J to be a real constant diagonal matrix with 
different eigenvalues. 

We will try to answer the question in the title first for the simplest class of Lax operators of 
Zakharov-Shabat type with real-valued J. In doing this we will be using the dressing method, one 
of the best known methods for constructing reflectionless potentials and soliton solutions. This 
paper is intended as a natural continuation of the work [3] published several years ago by two of 
the authors. 
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In Scction[3]below we first outline the well known facts about the soliton types of NLEE solvable 
by the sl(2) Zakharov-Shabat system. In Section [5] we treat the different one-soliton solutions 
for the sl(n) Zakharov-Shabat systems related to the subalgebras sl(p). Most of our results are 
illustrated for the si (5) system, but it is not difficult to extend them to any sl(n) system. The 
structure of the eigenfunctions of L(X) corresponding to the different types of solitons is outlined in 
Section[6] In the last Section we discuss possible generalizations to other Zakharov-Shabat systems 
having additional symmetry properties. The presence of symmetries modifies the spectrum of the 
scattering operator L. A typical example is the reduction of the type W(x, t, A*) = U(x, t, A). As a 
result the eigenvalues are pairwise symmetrical with respect to the real axis K, that is A + = (A - )* 
holds true. The soliton solutions are connected with two eigenvalues (doublet solitons). Another 
common situation is when we have a type reduction. The continuous spectrum of the Lax 
operator which is compatible with that reduction consists of 2h rays |13j . Thus the complex A- 
plane is split into identical sectors which possess equal number of eigenvalues. The soliton solution 
is associated with a multiplet of discrete eigenvalues (multiplet solitons). 

2. Preliminaries 

In this section we shall outline some basic features of the mathematical machinery we are about 
to use for the classification of soliton solutions. 

Integrability or more precisely S-integrability of a NLEE means that the NLEE can be presented 
as a zero curvature condition 

[L(A),Af(A)]=0, (1) 
of two first order linear matrix differential operators L(X) and M(A) of the form 

Lip(x, t, A) = id x ip + U(x, t, X)ip(x, t, A) = 0, (2) 
Mi/>(x, t, A) ee idti> + V(x, t, X)ip(x, t, A) = ip(x, t, A)C(A). (3) 

The potentials U(x, t, X) and V(x, t, X) are typically chosen as elements of some semismple Lie 
algebra g (the fundamental solutions ip(x, t, A) belong to the corresponding Lie group G). We shall 
mainly deal with the algebra sl(n). 

Remark 1. The compatibility condition ([I]) means that the Lax operators L and M possess the 
same eigenfunctions. The matrix C(A) depends on the definition of Jost solutions. 

The compatibility condition HJ which must hold true identically with respect to A takes the 
form 

id x V - id t U + [U(x, t, A), V(x, t, A)] = (4) 

and it is valid for any choice of C(A). For simplicity we shall resrict our considerations on scattering 
operators of Zakharov-Shabat type (GZS) 

L(X)ip(x, t, A) ee id x i> + {q(x, t) - XJ)tp(x, t, X) = 0. (5) 

The matrix J is a real traceless diagonal matrix, i.e. a real Cartan element of sl(n), while q(x,t) 
is a matrix with zero diagonal elements. Since J is a real matrix one can introduce an ordering 
of its elements J\ > Ji > ■ ■ ■ > J„. By carrying out a gauge transformation which commutes 
with J, we can always take q(x,t) to be of the form q{x,t) = [J,Q{x,t)\, i.e. qjj = 0. The linear 
subspace in sl(n) of matrix-valued functions q(x, t) = [J, Q(x, t)] are known in literature to be the 
co-adjoint orbit of g passing through J. The co-adjoint orbits can be supplied in a natural way 
with a non-degenerate symplectic structure which makes them natural choices for the phase spaces 
M.j and Hamiltonian structures of the corresponding NLEE. 

The class of NLEE related to L{X) are systems of equations for the functions Qjk(x,t), which 
may be written in the compact form 23, 20, 12,7: 



4 

id t Q + 2Y / A k [H k7 Q(x,t)} = 0, 
fc=i 



(6) 
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where H k , tr H k — are constant diagonal matrices and /(A) = X)fe=i ^ i s the dispersion law 
of the NLEE. Here and below we define: 

(adjX) ks = ([J,X}) ks = (J k - J s )X ks , ( a ,dj 1 X) ks = j^ T , (7) 

for all X £ A4j, i.e. Xkk = 0. The operator A is either one of the recursion operators A±, acting 
on the space M.j of n x n off-diagonal matrix- valued functions as follows 

A ± X = adf (id x X + P [q(x),X(x)}+iY^[Q(x),E k , k ] [ dy tr (Q(y), X{y)], E k , k ) ) . (8) 

V fe=l "'i 00 / 

where Po' is the projector adT 1 adj-. Choosing H\ = I = diag (Ii, . . . , I n ) ^ i, so that the 
dispersion law /(A) = —XI is a linear function of A we get a system, generalizing the well known 
N-wave equation: 

i[J,Q t ]-i[I,Q x }-[[J,Q],[I,Q}] = 0, (9) 

which contains N = n(n — 1) complex- valued functions Qij(x,t). 

In order to describe the soliton solutions we shall use the so-called dressing procedure [27] . For 
that purpose we need some basic facts on the direct scattering problem of L operator. 

Let ip±( x i t> ^) are two fundamental solutions of the GZS system ([5]). If they satisfy the require- 
ment 

lim ifj ± (x,t,\)e zXJx = i (10) 

they shall be called Jost solutions. The Jost solutions are interrelated via 

^(x,t,X) = iP + (x,t,X)T(t,X), (11) 

where T(t, A) is called a scattering matrix. The scattering matrix is x-indepedent and its time 
evolution is driven by the linear equation 

id t T+[f(X),T} = 0. (12) 

For the case of the iV-wave equation we have 

id t T-X[I,T(X,t)] =0. (13) 

Thus, if Q(x,t) satisfies the iV-wave system (j9|) we get: 

d t T kk {X) = 0, id t T jk {X) - X(I 3 - I k )T jk (X, t) = 0. (14) 

The set of matrix elements of T(A, t) must satisfy a number of relations. Indeed, they are uniquely 
determined by Q(x,t), i.e. by n(n — 1) complex functions of x, so it seems natural that there 
shouldn't be more that n(n — 1) independent functions among Tj k (X) for A on the real axis. Of 
course T(X,t) must satisfy the 'unitarity' condition detT(A, t) = 1. The rest of these relations 
follow from the analyticity properties of certain combinations of matrix elements of T(A, t). These 
analyticity properties must follow naturally from the corresponding fundamental analytic solutions 
(FAS) x ± ( x i t, A). 

The Jost solutions are well defined only for A 6 M, i.e. they do not have necessarily analytic 
properties beyond the real axis. This can be seen easily if one reformulates the problem ([5]) in 
terms of a Volterra type integral equation 

C ± (x, t,X) = l + i f dye iXJ (y- x U(y, m ± (y, t, X)e iXJ ^\ (15) 

J±oo 

where £,±(x,t, A) = ip±(x,t, X)e lXJx represents another set of fundamental solutions but this time 
to the linear problem 

id x £(x, t, A) + q(x, t)£(x, t, A) - X[J, £(z, t, A)] = 0. 

It is easy to see that only the first and the last columns of ifi+(x, t, A) and ip-(x, t, A) allow analytic 
extensions in A off the real axis; generally the other columns do not have analyticity properties. 
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Nevertheless it is again possible to introduce FAS [22"I[2"4| . Taking into account the ordering intro- 
duced above one is able to construct new fundamental solutions 



Ski + i dye^- J ^y- x Hq^) kl (y, t,X),k<l, 
if^dye^-^y-^iq^MvAX), k > I. 



(16) 



to possess analytic properties in the half planes C± of the spectral parameter. This definition can 
be rewritten using Gauss factors of the scattering matrix T 

X ± {x,t,X) = V-(x,t,A)S ± (t,A) = i> + (x,t,X)T*(t,X), (17) 

where T(t, A) = T T (t, A)(S ± (t, A)) -1 and x ± (x,t,X) = ^(x,t, X)e~ iXJx . The matrix elements of 
T ± (A) and S ± (A) can be expressed in terms of the minors of T(A). Here we note that their diagonal 
elements can be given by: 



mt_ 1 (A), 



T|(A)=m- -(A), 



TJ.(A)=mt(A), 
It. (A) = m'+^iX), 



(18) 
(19) 



where — = 1 and by mt(X) (resp. m k (A)) we have denoted the upper (resp. lower) 
principal minors of T(A) of order k, e.g.: 



m+(A) 



1 2 
1 2 



fc = l,...,5 



5-fc + l 5-fc 
5-fc + l 5-fc 



h h 



3k 



det 



T(A) 



T ■ T 



2 . 
2 . 

• 1}, 



T(X) 



I2jk 



(20) 
(21) 

(22) 



T ) 



As a consequence of the analyticity of the FAS, it follows that the minors rnt(X) (resp. ml" (A)) 
are analytic functions for A 6 C+ (resp. for A G C_). 

One can construct the kernel of the resolvent of L(X) in terms of the FAS [T2TI7] from which 
it follows that the resolvent has poles for all A^ which happen to be zeroes of any of the minors 
rrit(X). Therefore what we have now is that each of the minors rn^(A) may be considered to be 
an analog of the Evans function, and thus now, there is more than one Evans function. 

There exist different methods to solve a NLEE possessing a Lax representation : Gel'fand- 
Levitan-Marchcnko integral equation, Hirota method, dressing method etc. We shall use the 
dressing Zahkarov-Shabat method [27]. Let ipo(x,t,X) be a fundamental solution of Zakharov- 
Shabat's system with a known potential Uo(x,t, A) = qo(x,t) — A J. Consider a new function 
i>(x, t, A) = u(x, t, X)tpo(x, t, A) which is a solution to a Zakharov-Shabat's problem with some 
potential q(x, t) — A J to be found. This requires that u(x, t, A) satisfies 



qu — uqo — A [J, u] = 0. 



(23) 



The dressing procedure transforms the Jost solutions ^±,0(2^ tX), the scattering matrix To(t, A) and 
the fundamental solution % A) of the generalized Zakharov-Shabat system with a potential 
Uq(x, t, A) in the following fashion 



tp±(x,t,X) = u(x,t,X)iJj±fi(x 1 t,X)u ± 1 (X), 

T(t,X) = U+ (A)T (t,A)w: 1 (A), 
X ± (x,t,X) = u(x,t,X)xo(x,tX)uZ 1 (X). 



(24) 
(25) 
(26) 



The normalizing factors u±(X) — lim^^-too u(x, t, A) ensures the proper assymptotics of the dressed 
solutions ip±(x, t, A). 

Zakharov and Shabat t 24] proposed the following ansatz for the dressing factor u(x,t, A) 



u(x,t,X) = ± + (c(X)-l)P(x,t), c(A) = 



A- A+ 



(27) 
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where P is a projector (P 2 — P) which can be expressed via the fundamental analytic solutions 
(FAS) and A + (resp. A - ) is an arbitrary complex number in the upper (resp. lower) half plane C + 
(resp. C_). In the simplest case when rankP = 1 it reads 

\n(x, i))(m(x, t)\ 
P(x,t) = \\ . | / ' / ' 28 
v ' ; (m(x,t)\n(x,t)) K ' 

where 

\n{x,t)) =xt{x,t,\ + )\n ), {m(x,t)\=(m \{xo{x,t,\-))- 1 . (29) 

By taking the limit A — ► oo in equation (|23j) we obtain an interrelation between the seed solution 
(jo and the new one 

q = q + (\- -\ + )[J,P}. (30) 

Thus starting from a known solution of the NLEE we can find another solution by simply dressing 
it with some factor u(x,t, A). An important particular case is when go = 0- The dressed solution 
is called a 1-soliton solution. The fundamental analytic solution in the soliton case is given by a 
plane wave Xo( x ,t) = e ~ tX ( - Jx+It \ Repeating the same procedure one derives step by step the 
multisoliton solution of the corresponding NLEE, i.e. 

-> q ls -» q 2s -> . . . -» q ns . (31) 

Many integrable equations correspond to Lax operators that obey some additional symmetry 
conditions of algebraic nature. That is why it is worthwhile to outline some aspects of the theory 
of such Lax operators. 

Let an action of a discrete group Gr to be referred to as a reduction group be given on the set 
of fundamental solutions to the generalized Zakharov-Shabat system ([5]) as follows 

$(x,\) = Kip(x,k(X))K-\ (32) 

where k : C — » C is a conformal map. This action yields another action on the potential in the 
scattering operator L 

KU(x, /c(A))isr 1 = U(x, A) (33) 

A common case is when Gr = Z2 . Then the action of Z2 might involve external automorphisms 
of SL(n) as well 

ip(x, \) = K (ip T (x, fci(A)))" 1 K~ l , => KU T (x, hiX^R- 1 = -U(x, A), (34) 
${x, A) = Kij>*(x, k 2 {X))K-\ => KU*{x, k 2 {\))K- 1 = -U{x, A). (35) 

In particular, if Z2 acts trivially on the complex plane of the spectral parameter , i.e. k = id, then 
the symmetry condition (|34|) resricts the potential U(x, A) to a certain subalgebra of sl(n). For 
example, suppose K T = K then U(x, A) belongs to the orthogonal algebra so(n). The existence of 
a Z 2 reduction requires a modification of the dressing factor u(x, t, A) as follows 



u{x, t, A) = I + - I J P(x, t) + (c(A) - 1) P(x, t), (36) 

where P(x, t) is a projector of rank 1 and 

P(x, t) = KP T (x, t)K~ l . (37) 
The projector itself can be expressed through the FAS x ± ( a; i ^ ^) 

^'^tWtI' |n(a ! ,*))= X ( a: ,t,A-)|no), (m(x, t)\ = (m \( X (x, t, \ + ))-\ (38) 
{m(x, t)\n(x, t)) 
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3. Zakharov-Shabat system and sl(2) solitons 

The best known examples of NLEE are related to the Zakharov-Shabat system which is associated 
with the sl(2) algebra as follows 

Lip(x, t, A) = (id x + q(x, t) - \a 3 )tjj(x, t, A), q(x, t) = q + a + + q~a- = ( ®_ 9 Q J , (39) 

where a± — (cti ± ia-i)j2 and a\, 02 and 03 are the Pauli matrices. 

The class of NLEE for the functions q ± (x, t) related to (I39|) can be written in the compact form 

ia 3 d t q + 2f(A)q(x,t) = 0, (40) 

where /(A) is the dispersion law of the NLEE and A is one of the recursion operators, acting on 
the space A4q of off-diagonal matrix-valued functions as follows: 



4 2 



A ± X= l -[a^d x X}+ l -q{x) / dytr(q(y),[a 3 ,X(y)}). (41) 

J±oo 



The simplest nontrivial example of NLEE is related to a dispersion law of the type /(A) = — 2A 2 . 
This is the nonlinear Schrodinger equation 

iqt + qt x + 2(q+(x,t)) 2 q-{x,t) = 0, 

nT - - 2(q-(x,t)) 2 q+(x,t) = 0. (42) 
Another well known example is provided by a cubic dispersion law /(A) = 4A 3 , one gets the system 

4 + Qi xx + Qq + (t)q-(x, t)q+ = 0, 

It +lZx + ^(x,t)q + {x,t)q- =0. (43) 

directly linking to the Korteweg de Vries equation. 

As we discussed in the previous section the scattering theory is based on introducing Jost 
solutions of L(X), scattering matrix, fundamental solutions etc. In the sl(2) case the Jost solutions 
are 2x2 matrix-valued solutions defined by an analog of (I10p where the matrix J is simply 
substituted by 03. Then one introduces the scattering matrix T(\,t) by: 

T(V) = (V> + (M,A))-V-OM,A)= ( 6 +^ } "o-^f)' (44) 
which is a;-indepedent. The ^-dependence of the scattering matrix is driven by 

id t T+[f(\)a 3 ,T(\t)}=0, (45) 
Thus, if g ± (a;, t) satisfy the system of equations (|4"0")) we get 

d t a ± (X) = 0, t 6 ± (A) T 2/(A)& ± (A) = 0, (46) 

The matrix elements of T(A,i) are not independent. They satisfy the 'unitarity' condition 
detT(A) = a + a~ + b + b~ = 1. Besides the diagonal elements a + and a~ allow analytic exten- 
sion with respect to A in the upper and lower complex A-plane respectively. In fact the minimal 
set of scattering data which uniquely determines both the scattering matrix and the corresponding 
potential q(x) consists of two types of variables: i) the reflection coefficients (A) = 6 ± /a ± defined 
for real A £ K and ii) a discrete set of scattering data including the discrete eigenvalues Ajr € C± 
and the constants which determine the norm of the corresponding Jost solutions [19] . 

A simple analysis shows that the first column of ip+ allows analytic continuation in the lower 
half plane of the spectral parameter while the last one — in the upper half plane (for tp_ the 
opposite holds true) 

*/, + (x,t,\) = [1>l,tl>+], V-(M,A) = [^+ C]- (47) 
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The superscripts ± in the columns of the Jost solutions refer to their analyticity properties while 
the subscripts ± refer to different Jost solutions (with different limits of x). The fundamental 
analytic solutions are constructed in the following manner 

X + (x,t,\) = [Tp±,i>+], X~(je,t,X) = [tl>+,il>Z]. (48) 

The functions a (A) = detx ± (a;, A) are known as the Evans functions [26 , 3| of the system L(X). 
Their importance comes from the fact that they are ^-independent (see eq. (|46p L and therefore 
they (or rather lna ± ) can be viewed as generating functionals of the (local) integrals of motion. 
In addition it is known that their zeroes determine the discrete eigenvalues of L(X): 

a+(X+) = 0, A+GC+; a^A*) = 0, A.eC- (49) 

One can define the soliton solutions of the NLEE as the ones for which p ± (A) = for all AgK. 
Thus the soliton solutions of the NLEE (j40|) are parametrized by the discrete eigenvalues and the 
constants Cjjr whose i-dependence is determined from 

^=0, ^T2/X ± =0, g = f(Xt). (50) 

In fact we will analyze the various possible types of one-soliton solutions; in our case they are 
determined by one pair of discrete eigenvalues X^ £ C± and one pair of norming constants C^. 
Thus for the generic NLEE ([40|) we get just one type of one-soliton solutions. In order to derive 
its explicit form we shall use the dressing Zakharov-Shabat method [27] ■ In our case the dressing 
factor u(x,t, X) is given by a 2 x 2 matrix of the form (f2"7| where P is a projector of rank 1 (see 
formula ([28]) ) . Then the following relations hold 

P\n(x, t)) = \n(x, t)>, \n(x, t)) = ( ^ *j \ , (51) 

(m(x,t)\P(x,t) = (m(x,t)\, (m(x, t)\ = (m 1 (x,t),m 2 {x,t)) . (52) 
The transmission coefficients are transformed by the dressing procedure as follows 

a+(A)= C (A) a +(A), (A) = ^, (53) 

The sl(2) analog of formula ([30)) reads 

q(x, t) - q (x, t) = -(A+ - \-)[a 3 ,P(x, t)}. (54) 

By applying the above formulae to properly choosen constant vectors |no) and (mo| we can con- 
struct the eigenvectors of P(x,t) and as a result, obtain P(x,t) explicitly. It then remains only 
to insert it into eq. (p>4"]) in order to obtain the corresponding potential q(x,t) explicitly. It can be 
proved that the spectrum of L(X) will differ from the spectrum of Lq(X) only by an additional pair 
of discrete eigenvalues located at A ± € C± . 

A pure soliton solution is obtained by assuming qq(x, t) = 0; as a result we have: 

\n{x,t)) = e ( - l( * A++/+t)ff3 K), (m(x,t)\ = {m \e l{xX ~ + r t)a \ 
p (x,t) = — — — i„mx,t) „-* (x,t) ( 55 ) 



$ (x,t) = -i(X + - X~)x-i(f + -f-)t-1nK U 
$(x,t) = (X+ + X-)x+(f+ + f-)t, 

where and the constants K\ and K2 are given by: 



f ± = f(X ± ), K t = x m, i* = *H^. (56) 

Then the corresponding one-soliton solution takes the form: 

Q+(xt) - ^(A + -A-)e-^) (;i] _ «2(^-A-)e*M 

q [X ' >~ cosh<? (x,i) ' q [X,t) ~ K2 cosh<? (x,i) • [b ° 
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Remark 2. The eigenvalues A^ are two independent complex numbers, therefore in the denominator 
in eq. ([55]) we get cosh of complex argument. This function vanishes whenever its argument becomes 
equal to i(ir/2 + pir) for some integer p and the generic solitons of the (140|) may have singularities 
for finite x and t. 

One way to avoid these singularities is to impose on the Zakharov-Shabat system an involution, 
i.e. if we constrain the potential qo(x, t) by: 

q{x,t) = q f (x,t), i.e. q + = (q~)* = u(x, t). (58) 

Such constraint reduces the generic systems (|40|) to NLEE for the single function u(x,t); the 
second equation of the system becomes consequence of the first one. As a result eq. (|42|) becomes 
the NLS eq.: 

iu t + u xx + 2\u\ 2 u(x,t) = (59) 
while eq. (|43| goes into the MKdV-type equation: 

ut + u xxx + 6\u(x, t)\ 2 u x — 0. (60) 

This involution imposes constraints on all the scattering data; in particular we have: 

o+(A) = (o-(A*))*, b+(X) = (b-(X*)T, (61) 

^From the first relation above we find that the zeroes of the functions a (A) which are the eigen- 
values of Lq(X) must satisfy: 

A + = (A - )* = fi + iv, C + = (C-)*, P(x,t) = P\x,t). (62) 

So now the one-soliton solution corresponds to a pair of eigenvalues which must be mutually 
conjugated pairs. 

As a result we find that the expression for P(x,t) and the one for the one-soliton solution 
simplifies to 

1 / e #oo(x,t) e -i0 O i(x,t)\ 

P(X '* ) = 2 cosh *oo t) U* 0l(x ' 4) e-*o»(-.*) ) (63) 



$oq(x, t) = 2vx + 2ht - In 
$ 01 (x,t) = 2(ix + 2gt - argnj + argrip, 



'o 



where 



A* = n ± iv, / ± = g ± i/i. (64) 

Now both functions ^oo(^ji) and ^oiC^ji) become real valued. The denominator now becomes 
cosh of real argument, so this soliton solution is regular function for all x and t. 
One can impose on qa(x, t) a different involution 

q(x,t) = -q\x,t), i.e. q + = -(q~)* = u(x, t). (65) 

However it is well known that under this involution the Zakharov Shabat system L(X) becomes 
equivalent to an eigenvalue problem: 

£ip( x i t) A) = ia^dxip + a^q(x, t)ip{x, t, A) = Xip(x, t, A), (66) 

where the operator £ is a self-adjoint one, so its spectrum must be on the real A-axis. But the 
continuous spectrum of C fills up the whole real A-axis, which leaves no room for solitons. 
Finally, the Zakharov-Shabat system can be restricted by a third involution, e.g. 

q(x, t) = -q T {x, t), i.e. q + = -q~ = ~iw x . (67) 

Such involution is compatible only with those NLEE whose dispersion law is odd function /(A) = 
—/(—A). Therefore it can not be applied to the NLS eq.; applied to the MKdV eq. it gives: 

+ 6(w x (x,t)) w xx = 0, (68) 
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which can be integrated ones with the result v — w x : 

v t + v xxx + 6{v(x,t)) 2 v x — 0, (69) 

i.e. we get the MKdV eq. for the real- valued function v(x, t). It is well known also that the NLEE 
with dispersion law /(A) = (2A) _1 can be explicitly derived under this reduction and comes out 
to be the famous sine-Gordon eq. pQ: 

w xt + sm(2w{x,t)) = 0. (70) 



This second involution can be imposed together with the one in (|58f) . The restrictions that it 
imposes on the scattering data are as follows: 

o + (A) = (a-(A*))*, o+(A) = (o-(-A)), (71) 

Now if A + is an eigenvalue of L(X) then (A + )*, — A + and — (A + )* must also be eigenvalues. This 
means that we can have two configurations of eigenvalues: 

1. pairs of purely imaginary eigenvalues: 

\+ =iv= -(\+y, X" = -if = - (A - )*, (72) 

2. quadruplets of complex eigenvalues: 

A + = /i + iv — (A + )* = —fj, + iv, 

A - = /i — iv, — (A - )* = — fi — iv. (73) 

Thus we conclude, that the sine-Gordon and MKdV equations allow two types of solitons: type 
1 with purely imaginary pairs of eigenvalues and type 2 each corresponding to a quadruplet of 
eigenvalues. Type 1 solitons are known also as topological solitons, or kinks (for details see [4]). 
They are parametrized by two real parameters: v and |C + 1 so they have just one degree of freedom 
corresponding to the uniform motion. 

Type 2 solitons are known as the breathers and are parametrized by 4 real parameters: /i 
and v and the real and imaginary parts of C + . Therefore they have two degrees of freedom: one 
corresponds to the uniform motion and the second one describes the internal degree of freedom 
responsible for the 'breathing'. 

The purpose of presenting the above well-known facts in the above manner, was simply to make 
it clear that the structure, as well as the number of related parameters which determine what 
different types of solitons can exist, depend strongly on the type of, and the number of, different 
involutions that can be imposed on the system. 



4. TV- wave system related to si (3) 

In this subsection we are going to consider generic iV-wave system related to sl(3) along with 
its reductions. The 1-soliton solutions to this system are about to be derived as well. It proves 
to be convenient not to use a standard matrix notation but a notation which exploites the root 
structure of si(3), namely Qkn k,n = 1,2 stands for the component of Q associated with the root 
a = kct\ + ncti expanded over the the simple roots a,\ = e% — e2 and cti — e2 — e^. Taking into 
account that convention the s/(3)-iV-wave system consists of 6 equations of the form 

i(J% - J 2 )Qio,t - i(h - h)Qw,x + 3fcQiiQoY = 0, 
i(2J t + J 2 )Qu,t ~ i(2/i + h)Qu, x - 3kQ w Qoi = 0, (74) 
i(2J 2 + Ji)Qoi,t - i(2I 2 + h)Qoi,x + ZkQjnQn = 0, 

where k = J\Ii — I\Ji is an arbitrary constant. The rest of the 6- wave equations can be obtained by 
using the following tranformation: Qkn Qfc^> where Q-jz^ = Q—k,—n- This system can be solved 
via a dressing procedure with the dressing factor ((27)) . The 1-soliton solution obtained that way is 
given by the following expressions 
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Qu(z) = ^_^ e -^ + -+^-+-)>nX, (75) 
[m\n) 

Qoi(z) = ^—^e-^ + ^-^+^nlml 
(mm) 



where 



(m\n) = J^e~ l(A ~ A )zj njm^, z a = J a x + I a t, a = 1,2. 

3=1 

The other three fields can be derived from these by executing the following change of variables 

Qkn <-> Q^, e iX+Z] <-> e _lA ~% n 3 Q ^m 3 Q . 

Impose a Z2 reduction of the type 

KxU\x, X*)K{ 1 = E7(sr, A), =*> K X J*K^ = J, JCiQtJCf 1 = -Q, (76) 

where ifi = diag (ei, £2, £3) is an element of the Cartan subgroup H c SL{3) which represents an 
action of Z2 . This results in reducing the number of independent fields since we have 

Qiq = — ei£2<3*oj Qn" — — e i e 3<9ii, OoT = ~ e 2£3Qoi 
and therefore the number of equations from 6 to 3 as follows 

i(Jx - J 2 )Qw,t - i{h - h)Qw,x - 3fc£ 2 £3QnQoi = °' 
i(2Ji + J 2 )Qu, t - i(2h + h)Qn, w - 3fcQi Qoi = 0, (77) 
i(2 J 2 + Ji)Qoi,t - «(2/ 2 + /i)Qoi,3! - 3fc£i£ 2 Qi Qii = 0. 

The discrete eigenvalues of Z2-reduced operator L are complex conjugated, i.e. A~ = (A + )* = fi—if 
and the polarization vectors are interrelated via |n) = K\\m)* . The 1-soiton solution in this case 
is 

0j 1 , p v(z 1 +z 2 ) 
9,',/t>— Ml 

Om(V) - r -»M(2z 2 +zi) 2 / 3\* 

(nl^iln) 6 'W"o) • 

Remark 3. In general, the denominator of the expressions for the 1-soliton solution (|75j) can possess 
zeros for some x and t, i.e. we have singular solutions (exploding solitons). By imposing a certain 
reduction this effect can be annihilated. As it is seen from the solutions to the Z2-reduced problem 
with Ki = i we obtain a sum real exponents multiplied by some positive factors which do not 
vanish on the real axis. 

By imposing another Z2 reduction on the potential U(x,X), namely 

K 2 U T (x, -X)K 2 ; 1 = -U(x, A), K 2 J T K^ = J, K 2 Q T K^ 1 = Q. (79) 

where K 2 € H satisfies [Ki,K 2 ] = we obtain a Z2 x Z2-reduced si (3) iV-wave system. As a 
consequence we have a pair of purely imaginery eigenvalues X^ = ±if. Choosing K\ = K 2 = 1 we 
see that the three independent fields Qio(x,t), Qoi(x,t) and Qn(x,t) are purely imaginary while 
the polarization vector is real, |n)* = After introducing new variables 

Qio(x,t) = iq 10 (x,t), Q 01 (x,t) = iq x(x,t), Q n (x,t) = iq n (x,t), 

we derive a real 3- wave system for 3 real valued fields 

(Ji - J 2 )qio,t - (h - h)qio, x + 3fcqnq i = 0, 
(2Ji + J 2 )qu,t - (2/i + / 2 )qii,x - 3fcqi qoi = 0, (80) 
(2 J 2 + Ji)q i,t - (2/ 2 + ii)qoi,a! + 3fcqi qn = 0. 



How many types of soliton solutions do we know? 
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Since the dressing factor must satisfy the conditions 

(ut(a;,A*)) _1 = u(x,X), (81) 

(u T 0,-A)) _1 = u(x,X), (82) 

the projector P is real valued. In this case the discrete eigenvalues are purely imaginary, i.e. 
A* = ±iv. The 1-soliton solution is 

qll(x) = -2uP kl (x), P=W^, k?l 

(n\n) 

Taking into account that \n) = e" Jx \no) we derive explicitly the following result 





(nl) 2 - 


f- e 2uZ2 (nl) 2 + er 


-2v(z\+zi) 


(nl) 2 




K) 2 - 




-2v{z 1 +z 2 ) 


HY 


p2vzi 


K) 2 - 






(nl) 2 



qoiOM) : 

In the Z2 x Z2 case there exists another type of soliton solutions — these obtained by using a 
dressing factor of the form 



u(x, A) = 1 + (c(A) - l)P(x) + - K 2 P T {x)K^ 



= l + ^ P{x) + ^ KlP r (x)K ^ P W = MM. (84) 

These solutions are associated with 4 discrete eigenvalues of the scattering operator L: ±\ ± . In 
this sense they may be called quadruplet solitons unlike the solutions ([83)) which being associated 
with 2 eigenvalues ±iv represent doublet solitons. The vectors \n(x)} an d \m(x)) depend on the 
fundamental analytic solutions xti x i A) in the same manner as it is shown in (|29p . The dressing 
factor (|8~4")) is Z2 x Z2 invariant if the conditions hold true 

K x (u\x,\*)y X K^ 1 =u(x,\), (85) 
K 2 {u T {x,-\)Y 1 K^ 1 ^u{x,\), (86) 

provided that A + = (A - )* = u + iv and KxP^K^ 1 = P are satisfied. Moreover, we assume that 
the projectors P and K 2 P T K 2 1 are pairwise orthogonal, i.e. PK 2 P T = is valid. The resrictions 
on the projectors give rise to some algebraic relations on the polarization vectors, namely 

\m ) = K 1 \n }*, (n \K 2 \n ) = 0. 

To find the 1-soliton solution we take the limit A — > 00 in equation (|23|) and put qo = 0. Thus 
we obtain the following formula 

Q)% = (A- - A+)(P + K 2 P T K^) jk , j ± k. (87) 

Let K\ — K 2 — 1. Then Q* = —Q and using the above notation we have for the 1-soliton solution 

,1^21 



qio 



.44^ e K«i+») cosHzi _ Z2) _ 5l + 52]) 

(n|n)* 



= _^Kr|i e -^ cos[/i( 2 Zl + 22 ) - ^ + 5 3 ], 

(n|n)* 

= _4^jngf|| e _^ cog[ ^ (2z2 + _ ^ + 
{n\n)* 

\n)* = e 2 ^\nl\ 2 + e 2 » Z2 \n 2 \ 2 + e-^+^lnll 3 , n{ = |nj|e^. 
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5. Generalized Zakharov-Shabat system and sl(n) solitons 

For the sake of simplicity and clarity below, most of our discussions will be restricted to the 
case n = 5; however they also could easily be reformulated for any other chosen value of n. The 
corresponding Lax operator L(X) which is a particular case of eq. ([2]) with 

L = id x + U(x, t, A) = id x + [J, Q(x, t)] - A J, (89) 

/ Q12 Q13 Qu Q15 \ 
Q21 Q23 Q24 Q25 

Q31 Q32 g 34 Q35 , (90) 
Qn Q42 Q43 Q45 

V Q51 Qb2 Q53 Q54 / 



J = diag(Ji, J 2 , J 3 , J 4 , J5), Q(x,t) 



Furthermore, for definiteness we will assume that 

tr J = 0, Ji > J 2 > J 3 > 0, > J 4 > J 5 . (91) 

The M-operator in the Lax representation for the TV-wave equation ([9]) is given by: 

Mip(x, t, A) = id t ip + {[I, Q(x, t)} - \I)Tp{x, t, A) = -Xtp(x, t, X)I, (92) 

where / = diag . . . , I5) is a traceless matrix. 

As we discussed in the section Preliminaries the 1-soliton solution can be derived by using 
formula ([30]) 

q{x)= lim \(J -u{x,X)Ju- x {x,X)) = -(A+ - A")[J,P(x)]. (93) 

X— >oo 

where the projector P is of the form 

P[X)= M^Kff) ' K a; ))=^ + (^A+)K), (m(x)|=(m |xo(^A-). (94) 

The polarization vectors |no) and (mo| are constant 5-vectors. The 1-soliton solution is parametrized 
by: 

1. the discrete eigenvalues A* = // ± iv ; /i 1 * 1 determine the soliton velocity, determine the 
amplitude. 

2. the 'polarization' vectors. |no), (mo| parametrize the internal degrees of freedom of the soliton. 
Note that P(x) is invariant under the scaling of each of these vectors. Generically each 'polariza- 
tion' has 5 components, one of which can be fixed, say to 1. So each 'polarization' is determined 
by 4 independent complex parameters. 

We have several options that will lead to different types of solitons: 

- 1) generic case when all components of |no) are non- vanishing; 

- 2) several special subcases when one (or several) of these components vanish. The correspond- 
ing solitons will have different structures and properties. 

For the generic choice of |no) one finds: 

lim P{x,t)=P ± , P+=E n , P-=E nni (95) 

x — >±oo 

where the matrix E^j has only one non-vanishing matrix element equal to 1 at position fc,j, i.e. 

(Ekj)mp = fikmSjp- Therefore both the limiting values u±(X) and their inverse u±(X) are diagonal 
matrices: 

u+(A) = diag (c(A), 1,1,..., 1), u_(A) = diag(l,l,...,l,c(A)). (96) 
^From eqs. (|25|) for n — 5 we have 

Tij (A) = c(A)(To)y(A), j = 1,2,3,4; 

T JB (A) = (T ) j5 (A)/c(A), j = 2,3,4,5; (97) 

T i:j (X) = (T )„ (A), for all other values of (98) 



How many types of soliton solutions do we know? 
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This relation allows us to derive the interrelations between the Gauss factors of To (A) and T(A). 
In particular we find for the principal minors of T(X): 

m+(A) = c(A)m+ fc (A), (A) = m ok (X)/c(X), (99) 

where m k (X) (resp. TnZ(X)) are the upper (resp. lower) principal minors of T(X). Since \q {%■> ^ A) 
are regular solutions of the RHP then m t fe (A) have no zeroes at all, but eq. (fM]) means all mjf(A) 
have a simple zero at A = X^ . 

The generic one-soliton solution then is obtained by taking (x, i, A) = e — t H Jx + It ) As a result 
we get: 

(P(x,t)) ks = -i—ngm-e-*^'*-^*.), (100) 

n 

k{x,t) = 5^ng il mg )1 e- i ( A+ - x ~5*»^ > (101) 
P =i 

z k (x,t) = Jkx + ht, q% = -(\+-\-)(P(x,t))ks, (102) 

i.e. in all channels we have non-trivial waves. The number of internal degrees of freedom is 2(n — 
1) = 8. Note that the denominator k(x,t) is a linear combination of exponentials with complex 
arguments, so it could vanish for certain values of x, t. Thus the generic soliton (llOOp in this case 
is a singular solution. 

Next we impose on U(x, t, A) the involution: 

KU\x,t,X*)K~ 1 = U(x,t,X), if = diag(ei,...,e„), (103) 

with ej = ±1. More specifically this means that: 

Kq t (x 1 t)I<- 1 = q(x,t), Ku\x, t, X*)K~ 1 = u^ 1 (x, t, A), (104) 

and 

X+ = (X-)* =[i + iis, (mo | = (^K)) f . (105) 

Thus only |no) is independent. 

Then the one-soliton solution simplifies to: 

- - ^t^ ^^'''^"*' 1 ' (106) 

^rcd 

n 

fcredfot) - ^e p |<| 2 e 2t/ ^^. (107) 
P =i 

The number of internal degrees of freedom now is n — 1 = 4. If one or more of ej arc different, 
then this reduced soliton may still have singularities. The singularities are absent only if all ej are 
equal. 



5.1. Non-generic sl(2) solitons.. ^From now on we assume that the reduction (jl03|) with e p = 1 
holds. Here |no,i) has only two non-vanishing components. We consider here three examples with 
n = 5 and three different choices for the polarization vectors: 



1 



a ) K,i) 






\n 5 



b) |no) 



/ \ 

9 



'0 

V o / 



c) |n ) = 



/1\ 




V o y 



(108) 



In all these cases the corresponding one-soliton solutions q(x, t) are given by similar analytic ex- 
pressions, each having only two non-vanishing matrix elements: 



q k(x,t) = (q jk (x,t))* 
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iv(Jj — J fe ) e «( ar g( n o)^ ar g( n o)) e ^v(.-'j-.-'fc)(2:+M) J fct) 

cosh[i/(Jj — Jk)(x + Wjkt) + In \n 3 \ — In |n§|] 

where we remind that Wjk = {Ij — Ik)/(Jj — Jk), j < k. For the case a) we have j = 1, k = 5; in 
case b): j = 2, k = 4 and in case c) j = 1 and k = 2. 

The sl(2) soliton is very much like the NLS soliton (apart from the f-dependence) ; the NLS 
soliton has only one internal degree of freedom. 

The different choices for the polarization vector result in different asymptotics for the projector 
Pi(x,t): 



a) lim P(x,t) = En, hm P(x,t) = E 55 , 

x — *oo x — ►— oo 

b) lim P(x,t) = E 22 , lim P(x, t) = E Ai , 

x — *-oo x — > — oo 

c) lim P(x,t) =En, lim P(x,t) = E 22 , (110) 



x — > — oo 



In case a) the results for the limits of P{x,t) and for u±(X) are the same as for the generic 
case, see eqs. (|95|) . (|96|) . As a consequence, such sl(2) solitons requires the vanishing of all Evans 



functions m^(X) for A = X ± , see eq. (IM|) . 

In case b) from eq. (|25|) and from the appendix we get that such sl(2) soliton provides for the 
vanishing of m^(A) and m^(A): 

m 2 ( A ) = c ( x ) m o,i( X )> m tW = c(A)m^ 3 (A), 

m 2 (X) = mo i2 (\)/c(\), m 3 -(A)=m 0i3 (A)/ C (A), (111) 

whereas m^(A) = itiq 1 (X) and to^(A) = to^ 4 (A) remain regular and do not have zeros at A = A*. 
Likewise in case c) we get that only m^(A) and 7714 (A) acquire zeroes: 

m+(A) = c(A)m+ 1 (A), m 4 (A) = m^ 4 (A)/c(A), (112) 

and all the other Evans functions m^~(A) with j = 2, 3, 4, and m~ (A) with p = 1, 2, 3 do not have 
zeroes. 



5.2. Non-generic sl(3)-solitons. Here |no) has three non-vanishing components. We consider three 
examples of such polarization vectors: 



4 



a) |n ) = 





Wo/ 



b) 



( ° \ 

9 



\0j 



9 



"0 





V y 



Therefore the sZ(3)-solitons have two internal degrees of freedom. 
The asymptotics of the projector P(x, t) read as follows: 



a) lim P(x, t) = En, 

x — >oo 

b) lim P{x,t) = E 22 , 

c) lim P(x,t) = E n , 



lim P(x,t) = E 55 , 

x — > — OO 

lim P(X)t) = 

00 

lim P(x,t) = E 33 , 



X — > — OG 



a;— — 00 



(113) 



(114) 



Note that cases a) and b) in eq. (|114[) coincide with the corresponding cases in eq. (|108[) . 
Therefore the set of Evans functions that acquire zeroes will be the same as for the corresponding 
sl(2) solitons. In case c) of eq. (|114[) we have: 



m+(A) = c(A)m+ 1 (A), 
m 4 (A) = mn 4 (A)/cA), 



"4(A) = c(A)to+ 2 (A), 
m 3 (X) = m a3 (A)/c(A), 



(115) 
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whereas the remaining Evans functions m^~(A) with j = 3,4, and m~(X) with p = 1,2 remain 
regular. 

In case a) the corresponding one-soliton solutions acquire the form: 



a) q ls (x 1 t) = 



c) q ls (x,t) 



( q 13 q 15 \ 



9 13 935 



V915 935 / 

/ 912 923 0\ 

9i* 2 923 

9 13 92*3 



V 0/ 



b) q ls (x 1 t) = 



/0 0\ 

923 924 

92*3 934 

92*4 93*4 

\0 0/ 



where the matrix elements qk s (x,t) are given by: 
q ks {x,t) = (q sk (x,t))* 



| n l|2 e 2„(Ji2+/it) + | n 3|2 e 2 l /(J 3 x+Z 3 t) + | n 5|2 e 2KJ 5 



x+Int) 



and 



Jk = Jk - {Ji + h + Js)/3, h = h~ (h 



4)/3, 



Vks 



Jk + Js 



(116) 



(117) 



(118) 



This soliton has two internal degrees of freedom and is regular. 

Obviously it is by now clear how one can write down more complicated solitons like si (4) which 
would be characterized by polarization vectors of the form: 



a) K,i 



V o / 



b) K> 





W J 



(119) 



The si(4)-solitons will have three internal degrees of freedom. 

We note here that due to our choice of J in ((9i"|) . sZ(4)-solitons cannot give rise to generalized 
eigenfunctions. 



6. Eigenfunctions and eigensubspaces 



The structure of these eigensubspaces and the corresponding solitons becomes more complicated 
with the growth of n. 

In what follows we start with the generic case and split the 'polarization' vector into two parts: 



bo) = bo) + Mo); 



bo) 









nl 







n§ 


Mo) = 










nt 


W 




WoJ 



(120) 



and therefore 



b) = b) + M), b) =Xo f (^,^A+)b ), 



= X f (M,A + )Mo), 



(121) 



This splitting is compatible with eq. (j9T|) and has the advantage: if Xo { x ^i ^ + ) = e~' A+ Jx then 
b) increases exponentially for x — > oo and decreases exponentially for x — > — oo; \d) decreases 
exponentially for x — > oo and increases exponentially for x — > — oo, see also the lemma below. 
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What we will prove below is that one can take a special linear combination of the columns of 



Xq {x,t, A + ) which decreases exponentially for both 



and x — > — oo. Doing this we will use 

(122) 



the fact that 

X + (x,t,\+)\n ) = (l-P(x,t))x + (x,t,\+)\n ) = (i - P(x, t))\n(x, t)) = 0, 
Lemma 1. The eigenf unctions of L provided by: 

f+(x,t) = x + (x,t,X+)\ Po ) = - X + (x,t,\+)\d ), 
decrease exponentially for both x — > oo and x — > — oo. 



(123) 



Proof: ^From eq. (|122[) and (|120|) there follows that both expressions for f (x,t) coincide, so 
we can use each of them to ou r advantage, see eq. f|123|) . We will use also the fact that i — P{x, t) 
is a bounded function of both x and t. 

We start with 



lim f + (x,t) = lim X '- + (x,t,X+)\d ) = (1-P+) lim e- lX+iJx+I ^T-(X + )\d ), 



(124) 



where T~(A + ) is the lower triangular matrix introduced in eq. (|17p . If the potential is on finite 
support or is reflectionless then T~(A) is rational function well defined for A = A + . If the potential 
is generic then T~(A) does not allow analytic continuation off the real axis. Nevertheless T~(A + ) 
can be understood as lower triangular constant matrix (generalizing the constant Cq of the NLS 
case). Being lower triangular T~(A + ) maps \do) onto \d' ) = T^(A + )|do) which is again of the form 
(|120|) . i.e. its first three components vanish. Therefore 



lim e vax f + (x,t) = lim (1 - P+)e va 



I 



\ 



-i\+{JiX+I A t)±,l 

-i\+(J 5 x+I 5 t) 5,' J 



= o, 



(125) 



for any constant a > such that a + J± < 0. 
Likewize we can calculate the limit for x - 



-oo: 



lim f+(x,t) = - lim X' + (x,t,\l)\po) 

: — ' — oo x— > — oo 

= -(1-P+) lim e- lX+(Jx+It ^+(X + )\p ). 



(126) 



The upper triangular matrix § + (A + ) is treated analogously as T~(A + ). In the generic case it is 
just an upper triangular constant matrix which maps \po) onto \p' ) = § + (A + )|po) whose last two 
components vanish. Therefore: 



/ e -iA+(Jix+J 1 t) n l,/\ 



lim e" bx f + (x,t)= lim e" bx {l - P_) 



\ 



a -a+(j 2 x+j 2 t) n 2,/ 

-i\ + (J 3 x+I 3 t) n 3,/ 



J 



0, 



(127) 



for any constant b < such that J3 + b > 0. 
The lemma is proved. □ 

For the choices a) and b) of |no) in eq. ()108[) we define the square integrable discrete eigenfunc- 
tions using the splitting (|120p and eq. (|123[) . 



Remark 4- The choice c) for |roo,i) does not allow for the splitting (|120p . In this case we can 
introduce only generalized discrete eigenfunctions, / gcn (x,i), which are not square integrable. But 
upon multiplying by the exponential factor er vcx with c = (Ji + J2)/2, we can obtain square 
integrable functions f(x,t) = f (x,t)e~ ucx . See also the discussion in the next subsection. 
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The generalized eigenfunctions come up in situations when the splitting (|120|) is not possible, i.e. 
when either \p ) or \d ) vanish. Let us construct the generalized eigenfunction for the polarization 
vector |no) of case c) in eq. (|113p . Let (Ji + J 2 + J^l/i = a': then J[ = Ji — a', J 2 = Ji — a' and 
J'z = J3 — a ' are such that J[ > J 2 > J3 and J[ + J 2 + J 3 = 0. Let us assume for definiteness that 
J[ > J ' 2 > and > J3. Then we can split |no) into 



K> = \p' ) + K) 



\Po) 



( n \\ 




V / 



\d': 









v y 



(128) 



and define 

/ + ''(M) = X + M,A+)b ) = - X + (x,t,X+)\d' ), (129) 
Obviously f + ''(x,t) is an eigenfunction of the dressed operator L corresponding to the eigenvalue 
Then we can prove the following lemma: 



Lemma 2. The eigenfunction f +, '(x,i) is such that e Via x f '(x, t) decreases exponentially for 
both x — > ±00. 

Proof: The proof is similar to the one of lemma[T]and we omit it.D 

Since the polarization vector \uq) in case c) of eq. (|113|) does not allow the splitting (|120p 
the corresponding discrete eigenfunction will not be square integrable, so it will give rise to a 
generalized eigenfunction. 



7. Classification of solitons for a TV-waves system related to so(5) 

In this section we analyse how different kinds of reductions affect the classification of the soliton 
solutions to a nonlinear equation. This criterion is tightly connected with symmetries imposed on 
the auxiliary linear problem (the zero curvature condition). We shall consider in next subsections 
types of solitons which differ from one another in the number of eigenvalues associated with them: 
doublet solitons associated with 2 purely imaginery eigenvalues A = ±iv and quadruplet solitons 
associated with 4 eigenvalues situated symmetrically with respect to the real and the imaginery axis 
in C. This is the case when a Z2 x Z2 reduction is in action. Such type of reduction is compatible 
with the Lax representation of a NLEE to have a dispersion law obeying /(—A) = —/(A) (the 
TV-wave equation fulfills that restriction since /W-w(A) = —A/). 



7.1. N-wave system related to so(5). ^From now on we shall focus our attention on a TV-wave 
equation related to the so(5) algebra. This algebra has two simple roots Qi = e\ — e%, a 2 = e 2 , 
and two more positive roots: u\ + a 2 = e\ and a± + 2a 2 — e\ + e 2 = a max . When they come as 
indices, e.g. in Q a we will replace them by sequences of two integers: a — > kn if a = kct\ + na 2 . 
Moreover, we are going to use the auxiliary notation kn = —kati — na 2 . Thus the N-wave system 
itself consists of 8 equations. A half of them reads 

i{J\ - Js)Qio l t(x, t) - i(h - h)Qw,x{x, t) + kQu(x, t)Qj^{x, t) 
iJiQ\i,t{x, t) - ihQu,x(x, t) - fc(QioQ i + Qi2<2or)(£, *) 

i(Ji + J2)Qi2,t(x,t) - i(h + h)Qi2,x(x,t) - kQu(x,t)Q Q1 (x,t) 
iJ2Qoi,t( x i t) - ihQoi,x(x, t) + k(QjiQi2 + QwQn){x> t) 

The other 4 equations can be derived from those above by using the formal transformation Qkn *-* 
Qjn - C* ne ^ s arj ^ e t° integrate the system by applying the already discussed ideas — dressing method 



0. 
0. 
0. 
0. 



(130) 
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etc. For that purpose we make use of the dressing factor (j36|) . The 1-soliton solution reads 

A+-A- 



Qw(z) 



Quiz) = 



e i(X+ Z2 -X- Zl ) n l 2 



Iml + e^-^^ntml), 



A+ - A' 
(m|n) 



e -i\ z i n l TO 3 



Qa{z) 



A+-A- 



Qoi(z) = 



A+ - A" 



,-tA z 2 „2 3 



n m + e 



l o" l o ) > 



(131) 



fc=i 



i(A+-A )z k k k 
n rn Q ' 



jfe = 1,2. 



The other 4 field can be formally constructed by doing the following transformation 



Q 



kn 



Q 



kn 1 



Let us consider a 7L 2 reduction of the type KU'(X*)K _1 = £/(A) where iC = diag (ei, £2, 1, £2, ei) 
is an element of the Cartan subgroup (eft = ±1, k = 1,2). This means that Jk = J%, Q a must 
satisfy: 



Qto — — e i£2<2io> Qor ~ _£ 2<3oi, Qtt — — ei<3ii) Qti — ~ e i e 2<9i2- 

The corresponding NLEE is given by 4 equation 

i(J\ - J 2 )Qio, t {x,t) - i(h - I 2 )Qxo, x (x,t) - ke 2 Q n (x,t)Ql 1 (x,t) = 0, 
iJiQn,t{x,t) - ihQii, x (x,t) ~ k(Q 10 Q 01 + e 2 Qi2Q* m ){x,t) = 0, 
i(Ji + J 2 )Qi2,tOM) - i(h + h)Q\2,x{x,t) - kQu(x,t)Q i(x,t) = 0, 
iJ2Qoi,t(x,t) - ihQox,x{x,t) - kei(Q* n Qi2 + e 2 Ql Qii)(x,t) = 0. 

The Z2 reduction requires a dressing factor in the form 

u(x,t,X) = 1 + (^^y - 1^ P(x,t) + (c(X) - l)KSP*(x,t)(KS)-\ 

where P(x, t) is a projector of first rank (compare with (j3l)|) ) 

K\m*(x,t))(m(x,t)\ 



(132) 



(133) 



P(x,t) 



(m*[x,t)\K\m(x,t)) 



X- = (A+)* 



(134) 



(135) 



and S is the matrix of the metric in C which is involved in the definition of the orthogonal algebra 
so(5), namely 

so(5) = {A € sl(5) : A T S + SA = 0} , S l3 = (-ly^S^-j- 
The generic 1-soliton solution obtained by using the dressing method is the following 
2iv 



Qw(z) 



>.{z) 



ei 



(m*\K\m 

Qu(z) = 



^ 1 A*^,2 c i ( A+2 2-^ + *zi) 1 , f™^* 5 i(A+*z 2 -A+zi)\ 

{m ) m e v / + e 2{m ) m e v >\ 



2iv 



(ei(ml)*mle 



-i\ Zl 



- (™o)*™o e 



2iv 



(m*\K\m) 



m* \K\m 

( c ^ ls l*^, 4 - i ( A+ * 2 i+ A+z 2) , , ( rr ,2\* rn 5-i(\ + z 1 +\+*z 2 )\ 



Qoi(z) = 



2iv 



e2(m 2 )*m 3 e- lX+ * Z2 + (mg)*m^e^ A ' Z2 
3 1 2 



(m*\K\m) = ei\ml\ e 



l|2„-2^i 



(m*\K\m) 

e 2 \m 2 \ 2 e- 2 ^ + \m 3 \ 2 + e 2 \m*\ 2 e 2 » z * + e^m^e 2 ^ 



(136) 



where |mo) is a constant vector (polarization vector). By imposing certain resrictions on the com- 
ponents of \m(x, t)) we can obtain NLEE associated with some subalgebra of so(5). Let us consider 
several simple examlpes: 
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1. Suppose mj = TOq = 0. The only nonzero wave is Qoi(x,t) related to the simple root a 2 (of 
course, we mean an independent nonzero wave since Q- a2 is nonzero too). Thus we conclude 
that the solution is a sl(2) soliton. Another sl(2) soliton is derived when m\ = mjj = is 
satisfied. In this case Qn(x,t) is the nonvanishing component, repectively the sl(2) subalgebra 
is connected with the root e\ = a\ + a 2 . 

2. Let TOq = is fulfilled. Then we see that Qio(x,t) and Qi2(x,t) are nonzero waves. Since the 
corresponding Weyl generators commute this determines a representation of sl(2)(Bsl(2) w so(4) 
in so(5). 

3. Impose the resrictions (mj)* = wiq, (ra§)* = TOq and (mg)* = TOq. As a result we obtain 

Q w (z) = ^sinh26» cosh(^(zi +z 2 ))e" v(2:i ~ Z2+01 - 02) , 
^i 



(137) 



Qn(z) = - 2 ^\ W sinh 6> sinh(;/zi)e~' ( ^ 1+0l) , 
QiaC*) = ■^-sinh(20 o )cosh( I /(z 1 - za )) e -*/'(*i+* a +*i+* a ) j 

Q Q1 (z) = ^^ C osh6iocosh(^2)e^ (z2+02) , 
where we have used the representation 

3 3 

ml = ^Esmh0 o e #1 , to 2 , = cosh6» e^ 2 , O € K, 
v 2 v 2 

Z\i(a;, t) = 2 (sinh 2 9 sinh (vz-l) + cosh 2 9 Q cosh 2 (i/z 2 )) ■ 
In particular, if 9q = then one obtains a single wave 

Q 01 (M)= ^ e -^+^), ( 13 8) 
cosn(^z 2 j 

Remark 5. In the " soliton sector" the first two examples are trivial meaning that the 4- wave system 
(|133p is linearized. However, they have a nontrivial application when one constructs the 2-soliton 
soliton by the dressing 1 -soliton solution and when dressing a general FAS Xo( x i 



7.2. Doublet Solitons. In this subsection we are going to derive a 1-soliton solution to a 4- wave 
system with an additional Z 2 symmetry imposed on it. This is equivalent to a Z 2 x Z 2 symmetry 
condition imposed to it. Let the action of Z 2 x Z 2 in the space of fundamental solutions of the 
linear problem is given by 

X - (x, t, A) = K x (( X + )\x, t, A*)) _1 K- 1 

X~(x, t, A) = K 2 (( X + f(x, t, -A)) ^ K^ 1 

where S SO (5) and [Ki, K%] = 0. Consequently the potential U(x, t, A) satisfies the symmetry 
conditions ([76]) and (f79|) . The Z 2 x Z 2 -reduced 4- wave system reads 

(Ji ~ J2)qio,t(x,t) - (h - h)qio,x{x,t) + kqn{x,t)q 01 (x,t) = 0, 
Jiqi M (x,i) - hqxi >x (x,t) + k(q 12 (x,t) - q 10 (x,t))q Q1 (x,t) = 0, 

(Ji + J 2 )qi2,t(:M) - (h + h)qi2,x(B>t) - kq n (x,t)q Q i(x,t) = 0, 
J2<loi,t(x,t) - hqoi <x (x,t) + k(q 10 (x,t) + qi 2 {x, t))qu(x, t) = 0, 

where qio(x, t), qn(ic, t), qi 2 (a;, t) and qoi(£E, t) are real valued fields and their indices are associated 
with the basis of simple roots of B 2 introduced in the previous section, i. e. 

Qxo(x,t) = iq 10 (x,t), Qu(x,t) = iqn(x,t), Qiz{x, t) = iqi 2 (x, t), Q 0l (x,t) = iq i(x,t). 
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The constant k coincides with that one in the previous examples 

k := Jil 2 - J%I\. 

In accordance with what we said in previous chapter the dressing factor g(x, t, A) must be 
invariant under the action of Z2 x Z 2 , i.e. 

K x {ut{x,t,\*))~ l K t - 1 = u(x,t,X) (140) 

K 2 (u T (a;,t,— A)) -1 K^ 1 = u(x,t,X). (141) 

Let for the sake of simplicity require that K\ = K 2 = i. As a result we find that the poles of the 
dressing matrix are purely imaginery, i.e. 

A* = ±iu, v > 0. 

Thus the invariance conditions (|140p and (I141[) implies that the dressing matrix gets the form 

2iv 2iv 

u(x,t,X) = 1 + - P(x,t) - - -SP*(x,t)S, P*(x,t) = P(x,t). 

A — iv A + iv 

In the simplest case the explicit form of P[x, t) is 

= \m(x,t)){m(x,t)\ 
' (m(x, t)\m(x, t)) 

where the vector \m(x, t)) = e~ v ( Jx+It >\mo) is real. Therefore the solution turns into 
2v 



qw (x,t) = T 4 i T (e-^+^+^+^mlml + e^+AJx+fwy*]^^ 
(m\m) V / 

q n (x,f) = fe-^+^mJmg - e^+^mlml) , 

(m\ni) \ ' J 

qi2 (x,t) = ( e -Kii^*+(^/ 2 ) t ] m i m 4 + ^[(^^^+(fi-y^5 m 2j ^ 

q 01 (x,t) = (e-^+^mlml + e 

(771 m) V 



(142) 



m\m) \ 

These solutions can be rewritten in terms of hyperbolic functions as follows 

4v 

q 10 {x, t) = — — — N1.N2 cosh{zy[(Ji + J 2 )x + (h + I 2 )t] + S x + 8 2 }, 
(m\m) 

qn(M) = - 4t/ Niml amh[v(Jix + hi) +S{\, 
\m\m) 

Av 

qia(a?, t) = 7— ;— riViiV 2 cosh{i/[(Ji - J 2 )x + (h - I 2 )t] +5 X - S 2 }, 
(m\m) 

qoi(x, t) = . ^ . N 2 ml cosh[v(J 2 x + I 2 t) + <5 2 ], 
(m I to) 

(m(x,t)|m(s,t)) = 27V 1 2 cosh2( ; y(Jia; + /i7;) + <5i) + 2N% cosh2(v(J 2 x + I 2 t) + S 2 ) + (ml) 2 , 
where we have implied that m§ > for k — 1, 2, 4, 5 and therefore the following expressions 

5 1 ™4 



(143) 



5i:=-ln— 5 2 :=-ln — §, iVi = y mJ m 0! ^2 = V TO o TO o- 
^ tt7q 



make some sense 

'"' "~*Q = TOq ailU 7/Iq = TOq 



In particular, when mj = TOq and ttiq = TOq or in other words Si = 5 2 = we obtain 



qio(x,t) = -t— sinh(26» ) coshz/[(Ji + J 2 )x + (7i + J 2 )f] 
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qn(a;, t) = — ^— sinh#o smhu(JiX + lit), 
qi 2 (x,t) = -^sinh(26» )cosh^[(Ji - J 2 )x + {h - I 2 )t], 

qoi(^ji) = ~~r~~ cosh 9q cosh v(J 2 x + I 2 t). 
Ad 

where 9q £ R and 

Ad(x, i) := 2 (sinh 2 6» sinh 2 i/(Jia; + /it) + cosh 2 do cosh 2 i/( J 2 x + I 2 t)) . 

7.3. Quadruplet Solitons. This time to ensure the proper Z 2 x Z2-invariance of u(x, t, A) we add in 
the expression for it two terms more. The requirements (|140j) - (|14ip lead to the following dressing 
matrix 

u{x,t,\)-l + — —+ x _ ix+ y aTa^ , 

v ; (144) 

K 1 K 2 A*(x 1 t)(K 1 K 2 )^ 1 
A+(A+)* 

Taking into account the explicit formula (|144|) one can derive in the soliton case Qo(x,t) = 
the following relation 

[J,Q](x,t) = [J, A + KiSA*SKi - K 2 SASK 2 - K 1 K 2 A*K 2 Ki](x,t). (145) 

Like in previous considerations we decompose the matrix A(x, t) using two matrix factors |X (a:, t)) 
and \F(x,t)} and derive some differential equation for \F(x,t)) which leads to 

\F(x,t)) =e tX+iJx+It) \F Q ). 

The linear system for \X(x,t)) in this case is following 

1 K 2 SA(x,t)SK 2 K 1 SA*(x,t)SK 1 K x K 2 A*{x, t)K 2 K x \ 
1 2A+ + 2^ 2^ ) t}> = °- 

Starting from the above mentioned equation we can derive the following auxilliary linear system 

S|ffKtlH 2ffi Kt m in+ («) w 

2fi 2iv 2(A+j* 

where we introduced auxiliary entities 

\Y(x,t)) :=K 2 S\X(x,t)), \Z(x,t)) := KiS\X*)(x,t), \W(x,t)) := KiK 2 \X*{x,t)), 

\G(x,t)) :=K 2 S\F(x,t)), \H(x,t)) := A^|F*0M)>, \N(x,t)) := K x K 2 \F*{x,t)). 
To calculate \X{x,t)) we just have to solve the linear system shown above. The answer reads 



(\X)\ / a* b -c\ /S\F)\ 

\ Y ) \( xt )- a * ° ~ C 6 SlG) 

\Z) [x - t >~ A(x,t) \-b-c o S\H) 

\\W)J \-c-b a / \S\N)J 



(x,t), 
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where 



A(x, t) := \a(x, t)\ 2 + b 2 (x, t) - c z (x, t), a(x, t) 



(F(x,t)\S\G(x,t)) 



c(x,t) 



2A+ 

(F(x,t)\S\N(x,t)) 



liv 2fi 

Finally putting the result for \X(x,t)) in Q145P we obtain the generic quadruplet solution to the 
4- wave system associated with the B2 algebra (suppose K\ = K 2 = 1) 



qio(M) 
qn{x,t) 

qi2(a;,t) 

qoi(x,t) 
where 
a(x, t) 



Im 
Ini 
Im 
Im 



imN? 

a*Nx coshtipx + ip 2 ) (//cosh^T + tp 2 ) — ivcoshftp^ — ip 2 )) 

a*N x sinh(^i) - N{ sinh(y>*) 



a*N\ cosh((y9i — ip 2 ) (/icosh(iy9* — tp 2 ) — ii/cosh^* + tp 2 )) 

a*N 2 cosh(ip 2 ) — ■ N 2 cosh^j) 



No 



No 



1 



H -\- iv 



F 2 

Ni cosh 2ip! + N 2 cosh 2ip 2 + 



, , x mix, t) , , m(cc, t) 

b(x, t) = K . ' , c(x, t) = v ' ' , 



IV 



i(x,t) = |iVi| 2 cosh(2Re (fi) + \N 2 \ 2 cosh(2Re tp 2 ) 



,3|2 



N a := x /mZm 6 - a , 



(p a (x, t) := i\+(J a x + I a t) + i log a =1,2. 

2 m Q 



8. Discussion and further studies 

Here we shall outline some further topics which could be studied and which could lead to a deeper 
understanding of these soliton properties. 

The first obvious remark is that sl(n) contains as subalgebras also so(p) and sp(p) subalgebras. 
So it will be interesting to specify the conditions under which L(X) has solitons of type so(p) or 
sp{p). 

Second remark of the same nature is that one can start with L(X) related to so(n) or sp(n) 
algebras; such generalized Zakharov-Shabat systems allow one to solve special types of TV-wave 
systems whose soliton solutions have not yet been classified. Such systems, due to the additional 
symmetry, have a richer structure. 

The explicit form of the corresponding N-w&ve system related to these algebras has been re- 
ported in [5l l25llT5] . see also [9|fl8] . What could be done is to analyze the structure of its soliton 
solutions [9, 18J which are more involved due to the additional orthogonal symmetry involved. How- 
ever this symmetry complicates the construction of the dressing factors. Nevertheless, interesting 
new types of integrable cubic interactions could be obtained. 

Also, even more complicated types of solitons will be related to projectors of higher rank. The 
projectors P(x, i), which we used above, were all of rank 1. The rank 2 projector P 2 can be defined 
as: 



P 2 {x,t) = J2\ n k)M ks (x,t){n* s 



M sk (x,t) = (n*\n k ), 



M = M~ 



(146) 



Now each soliton will be parametrized by two polarization vectors; the corresponding eigensubspace 
will be two-dimensional too. Among the various types of rank-2 one-soliton solutions, there will be 
various possible configurations for the two polarization vectors. An example of a dressing factor 
u(x, t, A) constructed by a projector of second rank is the following one 



u(x, t, A) = i 



1 



-1 P 2 + (Vc^A) - l)P 2 
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Such type of a dressing factor was used by Wadati and coauthors [17] to derive the soliton solutions 
to a multicomponent Schrodinger equation relate to symplectic algebra sp(A). 

It is known in general how the machinery, well understood for the AKNS system such as Wron- 
skian relations, expansions over 'squared solutions', etc. can be generalized also for these types 
of systems. The dressing method, after some modifications, can also be applied, leading to the 
derivation of their soliton solutions. 

An interesting problem is the study of how the different possible reductions (see e.g. [9]) of these 
systems will influence the number of one-soliton types. 

Soliton interactions for various different types of solitons of these systems also present interesting 
problems. ^,From the results known for the iV-wave systems [23,20 it is known that new effects in 
soliton interaction, such as soliton decay and soliton fusion may arise. 
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